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Abstract. Development of fundamental physics relies on the constancy of various fundamental quantities such as 
the fine structure constant. Detecting or constraining the possible time variations of these fundamental physical 
quantities is an important step toward a complete understanding of basic physics. High quality absorption lines 
seen in the spectra of distant QSOs allow one to probe time variations of several of these quantities. Here we present 
the results from a detailed many-multiplet analysis, to detect the possible variation of fine-structure constant, 
performed using high signal-to-noise ratio, (~70 per pixel), high spectral resolution {R >45000) observations of 23 
Mg II systems detected toward 18 QSOs in the redshift range 0.4 < z < 2.3 obtained using UVES at the VLT. We 
validate our procedure and define the selection criteria that will avoid possible systematics using detail analysis of 
simulated data set. The spectra of Mg ll doublets and Fe ll multiplets are generated considering variations in a and 
specifications identical to that of our UVES spectra. We show our Voigt profile fitting code recovers the variation 
in a very accurately when we use single component systems and multiple component systems that are not heavily 
blended. Spurious detections are frequently seen when we use heavily blended systems or the systems with very 
weak lines. Thus we avoided heavily blended systems and the systems with Fe ll column density < 2 x lO^'^ cm"'^ 
while analysing the UVES data. To make the analysis transparent and accessible to the community for critical 
scrutiny all the steps involved in the analysis are presented in detail. The weighted mean value of the variation in 
Q obtained from our analysis over the redshift range 0.4 < z < 2.3 is Aa/a = (-0.06 ± 0.06) x lO"'^. The median 
redshift of our sample is 1.55 and corresponds to a look-back time of 9.7 Gyr in the most favored cosmological model 
today. The 3cr upper limit on the time variation of a is -2.5 x 10"" yr"^ < (Aa/aAt) < +1.2 x 10"^"^ yr"^ To 
our knowledge this is the strongest constraint from quasar absorption line studies till date. 

Key words. Quasars: absorption lines - Quasars: individual; HE 1341 - 1020, Q 0122 - 380, PKS 1448 - 232, 
PKS 0237 - 23, HE 0001-2340, Q 0109-3518, HE 2217-2818, Q 0329 - 385, HE 1347 - 2457, Q 0453 - 423, 
PKS 0329 - 255, Q 0002 - 422, HE 0151 - 4326, HE 2347 - 4326, HE 0940 - 1050, PKS 2126 - 158, Q 0420 - 388, 
HE 1158 - 1843 



1. Introduction 

Contemporary theories of fundamental interactions, such 
as SUSY GUT or Super-string theories, treating gravity 
and quantum mechanics in a consistent way, not only 
predict a dependence of fundamental physical constants 
with energy (which has been observed in high energy 
experiments) but allow for their cosmological time and 
space variations(Uzan, 2003). Detecting or constraining 

Send offprint requests to: R.Srianand 

* Based on observations collected at the European Southern 
Observatory (ESO), under the Large Programme "QSO 
Absorption Line Systems" ID No. 166.A-0106 with UVES 
on the 8.2m KUEYEN telescope operated at the Paranal 
Observatory, Chile 



the possible time variations of fundamental physical quan- 
tities is an important step toward a complete under- 
standing of basic physics. In the framework of standard 
Big-bang models one can probe the evolution of various 
physical quantities over the elapse time using measure- 
ments that are performed at different redshifts (z). As the 
energy of atomic transitions depend on the electromag- 
netic coupling constant a [its well measured laboratory 
value is eV?ic=l/137.03599958(52) (see Mohr & Taylor 
2000)], its possible time variation will be registered in 
the form of small shifts in the absorption line spectra 
seen toward high-z QSOs (Bahcall, Sargent & Schmidt 
1967). Initial attempts to measure the variation in a 
were based on alkali-doublets (Wolfe, Brown & Roberts 
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1976; Levshakov 1994; and Potekhin & Varshalovich 1994; 
Cowie & Songaila, 1995; Murphy et al. 2001). This method 
(AD method) uses the difference in the wavelengths of 
the doublets originating from the same ground state (i.e., 
— *^P3/2 and — >^Pi/2 transitions). The con- 
straints on the variation in a, Aa/a , is obtained by as- 
suming that the measured difference in the wavelength 
centroid of the doublets is proportional to to the low- 
est order. The best constraint obtained using this method 
is Aa/a = (-0.5 ± 1.3) x 10"^ (Murphy et al. 2001). 

Like absorption doublets one can also use the cen- 
tral wavelength of multiple emission lines originating from 
the same initial excited state. For example, Bahcall et 
al. (2003) use the nebular O iii emission lines at A5007 
and A4959 originating from "'^D2 excited level and derive 
Aa/a = (—2.0 ± 1.2) x 10~^ based on 73 quasar spectra 
over the redshift range 0.16<z<0.70. 

Studies based on molecular absorption lines seen in the 
radio/mm wavelength range are more sensitive than that 
based on optical/UV absorption lines. Till now good con- 
straints are available for two systems, Aa/a = (—0.10 ± 
0.22) X 10-5 at z = 0.2467 and Aa/a = (-0.08 ±0.27) x 
10~^ a.t z = 0.6847 if one assumes constant proton g-factor 
(Murphy et al. 2001a). It is also possible to have joint 
constraints on the fundamental constants using molecu- 
lar lines (see Chengalur & Kanekar, 2003). However, such 
studies at high z are not available due to the lack of molec- 
ular absorption systems at high z. 

The most suitable and accurate method for measur- 
ing Aa/a at high redshift is called the many- mult iplet 
method (MM method). This is a generalization of the AD 
method and was introduced by Dzuba et al. (1999). Unlike 
the AD method, this method uses absolute wavelength 
measurements of numerous absorption lines from differ- 
ent species. This method has been shown to provide an 
order of magnitude improvement in the measurement of 
Aa/a compared to AD method (Murphy et al. 2003 and 
references there in). However, in the bargain, this is the 
most vulnerable method against a number of systematics. 

Accurate laboratory wavelengths (up to an accuracy 
of few mA) are available for most of the important tran- 
sitions that are regularly detected in QSO spectra (see 
Table. 1). The sensitivity coefficients of the different line 
transitions from different multiplets to the variation in 
a were computed using many-body calculations taking 
into account the dominant relativistic effects (Dzuba et 
al. 1999 & 2002). If one uses the parametrization of these 
authors then the accuracy to probe the possible varia- 
tions in a depends very much on how well one measures 
the wavelength of the absorption lines at high z. 

In simple terms, MM method exploits the fact that the 
energy of different line transitions varies differently for a 
given change in a. For example, rest wavelengths of Mg il 
doublets and Mg l are fairly insensitive to small changes 
in a thereby providing good anchors for measuring the 
systemic redshift. Whereas the rest wavelengths of Fe II 
multiplets are very sensitive to small variations in a. Thus 
constraining relative shifts between an anchor and differ- 



ent Fe II lines allows to study the variations in a. However 
the accuracy depends on how well the absorption line pro- 
files are modeled. 

The absorption profiles are usually modeled using mul- 
tiple Voigt profiles that are defined by column density 
(TV), velocity dispersion(&) and redshift in addition to 
the rest- wavelength of the species. In general the Voigt 
profile decomposition is robust. However, the profile de- 
composition is not unique when the absorption lines are 
heavily blended and/or if the signal-to- noise ratio in the 
spectrum is not high enough. In addition, in real data, 
small relative shifts between absorption lines from differ- 
ent species can be introduced due to various systematic 
effects such as chemical and ionization inhomogeneities in 
the absorbing region, incorrect wavelength calibration, iso- 
topic abundances, and atmospheric dispersion effects etc. 
It is usually argued that these random systematic effects 
can be canceled by using large number of measurements. 
Indeed, applying the method to 143 systems over a large 
redshift range 0.2 < ^; < 3.7, Murphy et al (2003) mea- 
sured Aa/a = (-0.54 ±0.12) x 10"^ 

The detailed work based on MM method has been per- 
formed by only one group using simultaneous fitting of 
absorption line profiles of species covering a wide range of 
ionization states. Their important results need to be con- 
firmed independently using a data set that is optimized 
for the purpose. This forms the main motivation of this 
work. One of the main criticism made to the MM method 
is the apparent lack of transparency in the complex analy- 
sis (Bahcall et al. 2003). Here, we perform similar analysis 
to the one used by Webb and collaborators on a homoge- 
neous and very high quality data set. To make the whole 
process accessible for critical scrutiny, we describe in detail 
the procedure, and the individual fits. The main idea is 
that it is better to derive strong constraints on a small, al- 
though statistically significant, sample of well chosen sys- 
tems rather than loose constraints on a larger sample of 
strongly blended systems. 

The details of the observations and the quality of the 
data used in the analysis are discussed in Section 2. In 
Section 3, we discuss simulations that are performed to 
validate our Voigt profile fitting code in extracting the 
variations in a. We also use the simulation results to con- 
struct an optimum selection criteria that will ensure a best 
possible detection limit in Aa/a . In Section 4, we list 
all the Mg II, Fe ii systems that arc selected in our sam- 
ple. Comments on individual Voigt profile fitting are given 
in Section 5. Discussion and conclusions are presented in 
Section 6. 

2. Observations and data sample 

The data used in this study have been obtained using 
the Ultra-violet and Visible Echelle Spectrograph (UVES) 
mounted on the ESO KUEYEN 8.2 m telescope at the 
Paranal observatory for the ESO-VLT Large Programme 
"QSO absorption lines". This corresponds to a homoge- 
neous sample of 18 QSO lines of sight suitable for study- 
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Fig. 1. Properties of our sample: Panel(a) shows the distribution of median signal-to- noise (S/N) ratio per pixel 
computed over 30A intervals. Typical S/N = 60-80 per pixel is achieved over the wavelength range 4000-8000 A. 
Panel (b) shows the distribution of measured relative velocity shifts between Fe iiA2344 and Fe iiA2600 lines in our 
sample. As the q coefficients are similar for these transitions the velocity shift is expected to be distributed around zero 
in the case of good wavelength calibrations. Thus the plot provides an internal consistency check for the wavelength 
calibration. We notice that the accuracy in the wavelength calibration is of the order of 1-3 mA in the wavelength 
range of interest. Panel (c) shows the distribution spectral resolution measured from the calibration spectrum. Typical 
spectral resolution achieved is larger than 44,000. Continuous and dotted distributions in these panels are for the blue 
(4000-5500A) and red (5500-8500A) spectral range respectively. Panel (d) gives the histogram of the number of sight 
lines at a given redshift in which Mg ii and Fe ii absorption lines are in the observed wavelength range. All the 18 
sight lines in our sample cover the redshift range between 1 and 2 uniformly. 



ing various properties of the inter-galactic medium over 
a redshift range 1.7—3.2. All the objects were observed 
in good seeing conditions with 1 arcsec slit width. The 
data were reduced using the UVES pipehne, a set of pro- 
cedures implemented in a dedicated context of MIDAS, 
the ESO data reduction package. The main characteristics 
of the pipeline is to perform a precise inter-order back- 
ground subtraction for science frames and master flat- 
fields, and to allow for an optimal extraction of the ob- 



ject signal rejecting cosmic ray impacts and performing 
sky-subtraction at the same time. Usually a given wave- 
length range is covered by more than 6 exposures and 
removal of cosmic rays and odd flux fluctuations in a 
given exposure is robust. The reduction is checked step 
by step. Wavelengths are corrected to vacuum-heliocentric 
values and individual ID spectra are combined together. 
Air-vacuum conversions and heliocentric corrections were 
done using standard conversion equations (Elden (1966), 
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Fig. 2. Illustration of the methodology and reliability of our Voigt profile fitting code: The panels in the 
left hand-side show the velocity plot of the simulated Mg il and Fe ii absorption lines (points with error-bars) and the 
best fitted Voigt profiles (solid curve). The panel in the top right hand-side shows the as a function of Aa/a . The 
best fitted value and Itr error is shown in the panel. This is consistent with the input value of Aa/a — 0. Lower right 
panel shows the relationship between the input and recovered values of Aa/a . Each point is a result obtained from a 
single realization (there are 30 of them for a given input Aa/a ). Typical errors in these measurements are shown in 
the bottom right corner of the panel. The dark circles with error bars are the weighted mean value of Aa/a obtained 
from 30 realizations. 



StumpfF(1980)). The slit was always oriented along the 
parallactic angle and calibration exposures were taken be- 
fore or after the scientific exposures. Thus most of the im- 
portant concerns raised by Murphy et al. , are minimised in 
our observations. As we are using minimization method 
to detect sub pixel scale shifts the error spectrum is very 
crucial for our analysis. Great care was taken in computing 
the error spectrum while combining the individual expo- 
sures. Our final error is the quadratic sum of appropriately 
interpolated weighted mean errors and possible errors due 
to scatter in the individual flux measurements. 



In all the cases, spectra cover the wavelength range of 
3000-lOOOOA. A first guess continuum is fitted with an au- 
tomatic continuum fitting code. However for the system 
that are eventually used in the analysis, a local continuum 
normalization is performed using lower order polynomials. 
In all cases we have used the flux in the core of the sat- 
urated lines to estimate the background subtraction un- 
certainties. Corrections to this are applied whenever it is 
needed. We notice that continuum fitting is not unique 
in the Lya forest region and in the IR region where at- 
mospheric lines are plenty. Therefore we have avoided the 
absorption lines that are falling in the lower wavelength 
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side of Lya emission from the quasar and lines that could 
be blended with atmospheric lines in our analysis. 

The properties of the spectra in our sample are summa- 
rized in Fig. n Typical S/N~ 60 — 80 per pixel (typically 
O.O35A wide) is achieved in the whole wavelength range of 
interest (see left hand top panel in Fig.^. This is a factor 
2 improvement compared to the one used earlier for the 
analysis of Aa/a (see Fig. 2 of Murphy et al. 2003). The 
spectral resolution measured from the lamp spectrum is 
typically > 44, 000. 

We have investigated the accuracy of wavelength cal- 
ibration using the lamp spectra that are extracted in the 
same way as the object. We notice that individual emis- 
sion lines are well modeled by Gaussian. This confirms 
that the instrumental profile can be modeled with single 
Gaussian while fitting the Voigt profiles. The median dif- 
ference between the extracted and the actual wavelengths 
in most of the settings are close to zero. The root mean 
square of this deviation is always smaller than 3 mAin 
all the settings. We confirm the robustness of the wave- 
length calibration using other consistency checks. In par- 
ticular, the wavelengths of Fe 11 transitions strongly de- 
pend on a but the relative position of the lines (apart 
from Fe iiA1607) does not depend much on a. This means 
that, whatever the variation of a is, the relative velocity 
shifts between Fe 11 transitions should be consistent with 
zero in the case of good calibration. This is precisely what 
is observed (Panel b in Fig. Note in addition that the 
Fe II rest wavelengths are known to very high accuracy 
(~0.1mA). In principle similar exercise can be performed 
using Si IV and C IV absorption lines. However, unlike Fe 11 
lines, C iv wavelengths are known only to 2 mA accuracy. 
Even in this case the velocity shifts measured are consis- 
tent with zero within the allowed laboratory wavelength 
uncertainties (Petitjean & Aracil, 2003). When fitting the 
narrow doublets and multiplets from a given species in in- 
dividual settings, we have allowed for some adjustments 
in the instrumental resolution. This seems to be an impor- 
tant step as we realized that the seeing was much better 
than the slit width during most of the observing run. The 
metal line systems in the reduced spectra are identified 
using standard procedure. 

The next step is to validate our Voigt profile fitting 
codes and have some estimate of various possible system- 
atic effects. In order to do this wc first simulate the spec- 
tra with specifications similar to a typical UVES spectrum 
and perform the analysis on them. The details of the spec- 
tral simulations and the data analysis process is described 
in the following section. 

3. Chasing the systematics using simulated data 

As one is trying to measure Aa/a through Voigt profile 
decomposition of absorption lines, various systematic ef- 
fects can cause false alarm detection of Aa/a . The effects 
of random systematics can be averaged out by combining 
large number of measurements. However, the best strategy 
will be to have an a priori estimate of various systematic 
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effects so that a proper selection criteria can be applied 
while choosing the sample to achieve the best possible re- 
sults. In this section we try to understand various effects 
using simulated spectra. Such an exercise is also impor- 
tant to see how well our Voigt profile fitting code recovers 
the input value of Aa/a that is used while generating the 
spectrum. 

3.1. Generation of absorption lines 

In this section, we simulate data for different assumed val- 
ues of Aa/a and apply our procedure to recover the input 
value of Aa/a . We mainly concentrate on absorption lines 
from the Mg 11 doublet and Fe II multiplets. We assume 
the ratio, iV(Mg ii)/A^(Fe 11), to be constant in individual 
components. In order to mimic real data, the Voigt profile 
of the absorption line is convolved with an instrumental 
profile and appropriate noise is added to the spectrum. We 
model the instrumental profile as a single Gaussian with 
(7 = 1.69 pixels (with pixel size of 0.03 A) appropriate for 
A/AA = 42,000 as achieved in a typical UVES observa- 
tions (note however that the actual resolution achieved is 
somewhat higher in the real data). 

The noise in the spectrum is the combination of photon 
noise, readout noise, and the residuals from background 
subtraction. While the former dominates the error in the 
continuum the latter two will be the source of error in the 
core of saturated lines. We notice from our UVES data 
that the root mean square of residual at the bottom of 
highly saturated lines is ~0.4% of the unabsorbed contin- 
uum flux. The readout noise added to a typical normal- 
ized spectrum, mimics a Gaussian distribution with zero 
mean and a standard deviation of 4 x 10~^. The photon 
noise spectrum is obtain using Poisson statistic. To ob- 
tain a spectrum with a given signal-to- noise ratio, (S/N), 
we have generated Poisson random numbers {P[{S/N)'^]), 
with mean equals to (S/N)"^ and then assigned photon 
noise in the continuum of normalized spectrum, to be 
(P[(5/iV)2] - {S/N)'^)/{S/N)^, having Poisson distribu- 
tion with a = (S/N)~^. In the absorbed portion of the 
spectrum this value is scaled by square root of the inten- 
sity. In all our simulations we use the signal-to-noise in the 
continuum to be 70, corresponding to what is achieved in 
our data. 

In order to simulate a spectrum incorporating a vary- 
ing fine-structure constant (az), we use the analytic fitting 
function given by Dzuba et al. (2002), 

uj = ujo + qx. (1) 

Here wq and uj are, respectively, the vacuum wave number 
(in the units of cm~^) measured in the laboratory and in 
the absorption system at redshift z, x is a dimensionless 
number defined by x = (az/ao)^ ~ 1 where ao refers to 
the present value of the fine-structure constant and a^ its 
value at redshift z. The numerical value of the parameter 
q that are obtained using many body relativistic calcu- 
lations (see Dzuba et al. 1999, 2002) for different species 
are listed in the fifth column of the TableQ] Originaly uj is 



6 Hum Chand, Raghunathan Srianand, Patrick Petitjean & Bastien Aracil: fine-structure constant 



Table 1. Summary of atomic parameters that are used in our analysis 



Species 


Ao(A) 


Ref.(Ao) 


ujo{cm """) 


q 


id 


/ 


Ref.(/) 


Mg I 


2852.96310(8) 


L 


35051.277(1) 


+86(10) 


i 


1.830 


A 


Mg 11 


2796.3543(2) 


L 


35760.848(2) 


+211(10) 


a 


0.6123 


A 




2803.5315(2) 


L 


35669.298(2) 


+ 120(10) 


b 


0.3054 


A 


Al 11 


1670.7887(1) 


M 


59851.972(4) 


+270(30) 


j 


1.88 


B 


Si II 


1526.70709(2) 


M 


65500.4492(7) 


+50(30) 


k 


0.133 


A 




1808.01301(1) 


M 


55309.3365(4) 


+531(30) 


1 


0.00208 


A 


Fe II 


1608.45085(8) 


P 


62171.625(3) 


-1200(300) 


c 


0.0577 


C 




2344.2130(1) 


Q 


42658.2404(2) 


+1254(150) 


d 


0.114 


D 




2374.4604(1) 


Q 


42114.8329(2) 


+1640(150) 


e 


0.0313 


D 




2382.7642(1) 


Q 


41968.0642(2) 


+1498(150) 


f 


0.320 


D 




2586.6496(1) 


Q 


38660.0494(2) 


+1520(150) 


g 


0.0691 


D 




2600.1725(1) 


Q 


38458.9871(2) 


+1356(150) 


h 


0.239 


D 



Reference: 



oscillator strength:- (A) Morton D.C. (1991); (B) Prochaska et al.(2001) (C) Morton D. C(2003); 
(D) Bergeson et al. (1996); 

wavelength:- (L) Pickering et al. (1998); (M) Griesmann U. & Kling R.(2002); (P) Pickering et al. (2002) 
(Q) Nave et al. (1991) 



defined as uj = luq + qix -\- q2y with y = (az/cto) ~ 1 and x 
as defined above(Dzuba et al., 1999). When Aa/a is small 
q is basically q\ + 2^2 • The table also gives the laboratory 
values of wavelengths (with standard isotopic composi- 
tion for different elements) (Aq), wave-numbers (wg) E^nd 
oscillator strengths (/) used in this study. References to 
various parameters used are also provided in this table. 

3.2. Voigt profile fitting code 

To measure the relative shifts between an anchor and Fe ii 
lines, we fit all the lines simultaneously. In addition, we 
have modified our Voigt profile fitting code (Khare et al. 
1997) to take into account the dependence of rest wave- 
lengths of different species on a. We improved the con- 
volution algorithm using a 41 point Gaussian quadrature 
integration scheme, lines in high signal-to-noise data. The 
FWHM of the Gaussian kernel is considered to be X^/R. 
Here Ar is the median wavelength over the region of ab- 
sorption line being fitted and R is the spectral resolution. 
We have taken a ±3cr (cr = FWHM/2.3548) wavelength 
range around each pixel of absorption line for convolution. 
The Gaussian quadrature integration method is adopted 
because it provides the most accurate estimate of the area 
under the Gaussian kernel that is used for normalizing 
the truncated Gaussian probability function. Our code 
also takes into account the variation of spectral resolution 
across the spectrum. 

3.3. Measurement of Aa/a in a single component 
systems 

In order to find the best Aa/a consistent with our data, 
we fit the absorption lines using the modified rest wave- 
lengths (see equation (1)), varying Aa/a in the range from 



-5.0 X IQ-^ to 5.0 X IQ-s in steps of 0.1 x 10"^. For each 
value we obtain Xmin varying N, b and z. The value 
of Aa/a at which Xmin minimum is accepted as the 
best possible Aa/a value, provided the reduced of fit 
is ~ 1. Following the standard statistical procedure (Press 
et al.(2000) see pages p690-691) we assign Ic errors to the 
above defined best value of Aa/a by computing the re- 
quired change in Aa/a so that Ax^ = — Xmin = f • To 
be on the conservative side we have taken the maximum 
change in Aa/a around the best value of Aa/a that is 
required to produce Ax^ = 1 as our error on Aa/a . 

We test the method discussed above using a simulated 
single component system (see FiglSJ. The profiles of Mg ii 
and Fe ii absorption lines that are generated for Aa/a — 
are shown in the left side column of the figure. In the top 
right side panel we show the best fitted as a function 
of Aa/a . The location of the minimum and la error in 
Aa/a are shown using a point with an error bar. In or- 
der to investigate the systematics in detail, we have sim- 
ulated single component systems with column densities 
in the range 10^^ to 10^^ cm^^ (we consider the num- 
ber of systems with a given A^(Mg ii) is proportional to 
iV(Mg (Srianand & Khare 1994)), velocity disper- 

sion uniformly distributed in the range 3 to 10 km s~^and 
input Aa/a varying between — 10~^ and +10~^. In addi- 
tion to that we have also varied the seed for the random 
number generator that produces the error spectrum. We 
use our Voigt profile fitting code to recover Aa/a in each 
case. The recovered value of Aa/a in each realization is 
plotted (dots) as a function of the input Aa/a in the lower 
right panel in Fig. |21 Typical error in an individual mea- 
surements is also shown in the plot. There are about 30 
realizations for a given value of input Aa/a with a varied 
N, b and noise spectrum. The dark circles with error bars 
are the weighted mean value of the output Aa/a using 



Hum Chand, Raghunathan Srianand, Patrick Petitjean & Bastien Aracil: fine-structure constant 



7 




r Aa - Acx 1 (10-^) 

L out *-Y inJ ^ 



Fig. 5. Effect of blending in the determination of Aa/a : Absorption spectra of Mg ii and Fe ii were simulated 
for given N, b, and spectral resolution similar to that of our data, introducing spectral shifts corresponding to a given 
value of Aa/a . Top panels show the relationship between the input and derived value of Aa/a in the case of a single 
clean component (left-hand side) and a blend of two components (right-hand side). A typical absorption profile is also 
shown in these panels. Dots are the values from individual realizations; the points with the error bars are the weighted 
mean over 30 realizations. The lower panels give the distribution of the recovered Aa/a around the true one. Single 
(left) and blended (right) cases are considered respectively. Fits by Gaussian distributions are over-plotted. 



all the 30 measurements. The weighted mean recovered 
values of Aa/a are consistent with the input values. 

It is clear from Fig. |21 that our procedure works well 
for the simple single component systems. We notice that 
the distribution of differences between Aa/a output and 
input values is well fitted by a Gaussian (see Fig. with 
fj — 0.23 X 10^^. Based on central limit theorem and 
Gaussian distribution of the recovered values around the 
mean, we expect the accuracy of the estimated mean value 
to be ^ (t/\/N for N measurements. Strong blending and 
internal structure in the cloud will imply larger uncer- 
tainty. 



The improvement we expect from the enhanced S/N 
ratio in our sample compared to other studies is visualized 
in Fig. 13 Here, we show the results obtained for individ- 
ual realizations of single component systems with input 
Aa/a = for S/N = 30 and 70 in the top and bottom 
panels respectively. The shaded region in the figure shows 
the mean and a in the distribution of recovered values. 
The distribution of individual values has a of 0.409 x 10~^ 
and 0.239x10"^ for signal-to-noise ratios 30 and 70 re- 
spectively. It is clear that the improvement in the S/N 
ratio gives a factor 2 improved accuracy in the recovered 
value compared to what has been done previously. 
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Fig. 3. The effect of signal-to-noise ratio: The fig- 
ure shows the measured Aa/a and la error bar from 30 
realizations of absorption systems all simulated with an 
input value Aa/a ~ 0. The shaded region is the standard 
deviation (la) in the measured values around the mean 
value (shown by horizontal dashed line). The absorption 
systems are simulated using random realizations of iV, 6 
and noise spectrum. We obtained a = 0.409 x 10~^ and 
0.239x10"^ based on 30 realizations with signal-to-noise 
30 and 70 respectively, ft is apparent from the figure that 
better S/N will improve the measurements of Aa/a. 

In all the simulations discussed till now we have used 
only strong but unsaturated Fe ii lines (i.e iV(Fe ii)> 
10^'^ cm"'^). This ensures that the absorption profiles are 
well defined. However, when the lines are weak then the 
flux in the line is appreciably affected by Poisson noise. 
One can therefore infer that the scatter in the measure- 
ments will be larger if weak lines are used in the analysis. 
We demonstrate this using simulations of weak lines. For 
simplicity we have assumed input Aa/a =0, A^(Mg ii) = 
1.5 X 10^2 cm-2 and A^(Fe ii) = 5x10" cm'^. We have 
simulated 30 systems by just changing the error spectrum. 
The results are summarized in Fig. 01 In the left hand 
side panel we show the profiles for one realization. From 
this figure one can appreciate the distortion in the ab- 
sorption profile of Fe ii lines caused by Poisson fluctu- 
ations. The top panel in the right hand side shows the 
recovered values of Aa/a . As expected the error bars 
are larger (one a of the distribution around the mean is 
2.92x10"^). Typical errors are then a factor 10 higher 
than that we derive with strong single component sys- 
tems. Thus if one is going to use 1/cr^ weighting of dif- 
ferent measurements, only one single strong and narrow 
component can be as important as about 100 weak sys- 



Fig. 4. Simulations of weak lines: The left hand side panels 
show the velocity plot of the absorption lines simulated in 
one of the realizations. The best fitted Voigt profiles are 
overplotted. The right-hand side panels give the recovered 
value of Aa/a in individual simulations for two values of 
N(Fe ii). In all cases the input Aa/a is zero and only the 
error spectrum is changed. The shaded regions in these 
panels gives the standard deviation around the mean. The 
numerical values are also given in each panel. 

terns. When considering 30 realizations, we find a weighed 
mean of Aa/a = (-0.828 ± 0.508) x 10"^ that is a 1. Go- 
deviation with respect to the zero input value. The results 
shown in the bottom right panel in Fig.0]are for A^(Fe ii) 
= 10^2 cni-2 and iV(Mg ii) = 3x10^2 cni-2. The distribu- 
tion of individual measurements has a — 1.43 x 10^^ and 
the weighted mean value Aa/a = (0.076 ± 0.273) x lO"'''. 
Thus this study clearly demonstrates that if the lines are 
weak and the profiles are dominated by Poisson errors then 
a marginally significant falls alarm detection of Aa/a can 
be obtained even when one uses large number of systems. 
This can be avoided in the analysis of real data by care- 
fully estimating appropriate lower limit cutoff in column 
density for the sample. Based on our simulations and data 
(see below) we obtain this cutoff for Fe il column density 
to be of the order of 2 x 10^^ cm"^ for the median S/N 
ratio corresponding to our data. 

3.4. Aa/a measurement in strongly blended systems 

Next we consider the case of strongly blended two com- 
ponent systems. Here N and b for the two components 
are chosen in the same way as for the single component 
systems (see section 3.1). The separation between the two 
components is taken between 3 and 5 km s~^. This means 
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Fig. 6. Results for well separated blends: Left hand- 
side panels show the velocity plot and the best fitted Voigt 
profiles to Mg ii and Fe ii absorption lines from one of the 
realizations. In all cases considered here the velocity sepa- 
ration between the components are always more than the 
individual values of the b parameters. In the right hand- 
side panels we compare the results obtained for heavily 
blended systems (top panel), well separated blends (mid- 
dle panel) and single component systems (lower panel). 
The standard deviation measured around the mean in dif- 
ferent cases is also given in the corresponding panels. It 
is clear that systems with well separated sub-components 
are as good as that of single component systems. 



that the subcomponents are separated by less than the 
width of individual components. The relationship between 
the input and the recovered Aa/a is shown in top right- 
hand panel of Fig. |S1 This panel also shows a typical sim- 
ulated profile. There are 30 realizations for a given input 
Aa/a . Typical error in the recovered values is shown in 
the figure. It is apparent that individual measurements 
deviate much more from the actual value compared to the 
single component case. Unlike in the case of single com- 
ponent systems the deviation with respect to the actual 
value can not be approximated by a single Gaussian profile 
(see bottom right panel in Fig.EJ. There is an extra tail in 
the distribution on both sides. This clearly demonstrates 
that statistically significant deviations from the real value 
is possible for a non- negligible fraction of the systems. To 
avoid this, we will predominantly use systems that are 
not strongly blending (i.e, for which the sub-component 
separations are larger than the individual b parameters). 



3.5. Aa/a measurements from well separated blends 

The next step is to investigate the case of multicomponent 
systems where subcomponents are marginally blended. 
That is the cases where the separation between individual 
components is larger than the individual b values. Results 
are summarized in Fig. [S] The left-hand side panel shows 
velocity profiles for one of the realizations. In the right 
hand-side panels we compare results obtained for the case 
considered here (middle panel) with that obtained for sin- 
gle component systems (lower panel) and strongly blended 
systems (top panel) . The a of the distribution around the 
mean are 0.221 xlO~^ 0.261x10"^ and 0.519x10"^ for 
the single, weakly blended and strongly blended systems 
respectively. It is apparent that the mildly blended case 
gives very similar results as that obtained with single com- 
ponents. 

3.6. Conclusions from simulations 

In short, we can conclude from these simulations that: 

1. Best constraints on Aa/a are obtained either from sin- 
gle component systems or well resolved multiple com- 
ponent systems. 

2. Increasing the signal-to- noise ratio from S/N — 30 to 
70 increases the accuracy of Aa/a measurements by 
about a factor of two. 

3. It is better to avoid weak lines while extracting 
Aa/a as their profiles can be distorted by Poisson 
noise. Thus, weak lines in the low signal-to-noise 
data can result in false alarm detections of non-zero 
Aa/a values. 

4. There is a non-negligible probability to derive a statis- 
tically significant deviation from the actual value when 
one considers highly blended systems (i.e., systems 
where the component separations are smaller than the 
individual b values). Thus it is better to avoid complex 
blends in the analysis. 

In what follows we use these results as guidelines for choos- 
ing the systems in our sample for the final analysis. 

4. Description of the LIVES sample 

The Mg II systems with detectable Fe il absorption lines 
present in our data are listed in Table. [21 The table pro- 
vides the QSO name, the emission redshift (zcm), the red- 
shift of the absorption line systems (zabs), the species 
that are detected (notations are as in Table. No such 
Mg II systems are detected along 2 sight lines (toward 
HE 1158 - 1843 at Zom=2.449 and Q 0420 - 388 at Zcm= 
3.177). Last column of the table is meant for comments 
on the systems. 

In total, there are 50 systems with detected Mg ii and 
Fe II lines at a detection limit of ~1.5 x 10^^ cm~^ for 
the strongest Fe iiA2382 line. The corresponding limit for 
Fe IIA2374 is 1.5x10^^ cm^^. In the following we discrim- 
inate systems with A^(Fe ii) < 2 x 10^^ cm~^ as simula- 
tions have shown that spurious effects can affect results in 
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Table 2. List of Mg ii & Fe ii systems 



QSO 



^abs 



Transitions'' Comments 



HE1341-1020 2.135 



Q0122-380 



PKS1448-232 
PKS0237-23 



Q0109-3518 
HE2217-2818 



Q0329-385 

HE1347-2457 

Q0453-423 



PKS0329-255 
Q0002-422 



HE2347-4342 
HE0940-1050 



2.190 



2.220 
2.222 



HEOOOl-2340 2.263 



2.404 
2.414 

2.435 
2.611 
2.658 



2.703 
2.767 



HE0151-4326 2.789 



2.871 
3.084 



PKS2126-158 3.280 



0.8728 
1.2778 
1.9154 
0.8225 
0.8585 
1.2433 
1.9102 
1.5847 
1.1846 
1.3650 
1.6371 
1.6575 
1.6724 
0.4524 
0.9489 
1.5855 
1.6517 
2.1839 
1.1828 
1.3499 
0.7862 
0.9425 
1.6917 
1.6277 
1.5556 
0.7627 
1.4379 
1.4392 
1.5082 
0.7261 
0.9083 
1.0394 
1.1492 
1.6302 
1.8584 
2.3004 
0.9926 
0.8367 
1.5418 
1.9888 
2.1678 
2.3018 
0.6632 
1.7319 
1.7962 
1.0598 
1.7893 
1.9182 
0.6631 
2.0225 



abdefghi 

abcdefghi 

abcdefghi 

abdfghi 

abefgh 

abdefi 

abfh 

abfh 

abfgh 

abcdefghi 

abcdefghi 

abcdefghi 

abcdefghi 

abdefghi 

abfh 

abicdefgh 
abfi 

abdefghi 

abdfghi 

abcdefghi 

abfh 

abdfghi 

abcdefghi 

abcdefghi 

abcdfghi 

abefghi 

abdhi 

abcdefghi 

abdfh 

abdefghi 

abdefghi 

abh 

abdefghi 

abf 

abdfh 

abcdef 

abdfghi 

abdefghi 

abdefghi 

abdefghi 

abcdefghi 

abcdefi 

abdefghi 

abf 

abdfgh 

abdgh 

abdefghi 

adefghi 

abdefg 

abdefgh 



weak 
weak 
weak 



weak 
weak 

weak 
weak 

weak 

weak 
Lya 

weak 
+ 

weak 

+ 

weak 
Lya 

weak 



weak 
weak 
no anchor 

+ 

no anchor 
Lya 



case the lines are weak. Therefore systems with N{Fe ii) 
Jess than this limit (15 systems) are marked as "weak" 
in the table and are not considered in our main analysis. 
However, for completeness, we present the results obtained 
from these systems as well in the concluding section. 

Two systems at Zabs^ 1.0598 and Zabs = 1.9183 toward 
HE 0940+1050 could not be considered in our analysis be- 
cause there are no good anchor lines in these systems. In 
the former system only Mg iiA2796 is present in the spec- 
trum. In the case of the Zabs=l-9183 system. Si ii and Al ii 
are blended with Lya lines and Mg ii is heavily affected by 
atmospheric lines. In addition, at Zabs= 0.726176 toward 
Q 0453-423, Zabs= 0.8367 toward Q0002-422 and Zabs= 
0.6631 toward PKS 2126-158, the Fe ii hues are heavily 
blended with Lya systems and are therefore marked as 
"Lya'. Thus we are left with 29 systems that have strong 
Fe II lines, good anchor lines and arc not contaminated by 
Lya or atmospheric absorption. These systems form the 
basic dataset that we use for measurement of Aa/a . 

In summary, based on the results from the simulations, 
we apply the following selection criteria to derive reliable 
Aa/a : 



1. We consider only lines with similar ionization poten- 
tials (Mg II, Fe II, Si II and Al ll) as they are most 
likely to originate from similar regions in the cloud. 

2. We avoid absorption lines that are contaminated by 
atmospheric lines. 

3. We consider only systems that have iV(Fe ii)> 2 x 
10^^ cm~^ which ensures that all the standard Fe ii 
multiplets are detected at more than 5cr level. 

4. We demand that at least one of the anchor lines is not 
saturated so that the redshift measurement is robust. 

5. We also avoided sub-DLAs (i.e iV(H l)> 10^^ cm'^) 
as these systems may have ionization and chemical in- 
homogeneities. 

6. We do not consider strongly saturated systems with 
large velocity spread (complex blends); however in 
such systems whenever we find a well detached satel- 
lite components we include these components in the 
analysis. 

7. finally, based on the component structure resulting 
from the Voigt profile fits of systems that are not com- 
plex blends, we retain only systems for which the ma- 
jority of components are separated from its neighbor- 
ing components by more than the h parameters. 



''Transitions as identified in Tabled 

"weak" refers to systems with Ai'(Fe ll) < 2 x lO^^cm"^. 

"+" indicates systems that are not used in the 

analysis (see discussion). 
"Lya " refers to systems which are blended in the Lya forest, 
"no anchor" refers to systems without good anchor lines. 



Application of the above conditions resulted in 23 systems 
on which to perform the measurement (six single compo- 
nent systems, six well separated doubles, six systems with 
three components with at least one well detached from the 
rest, and five systems with more than 3 components) . 
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5. Notes on individual systems in our sample 

5.1. Systems along the line of sight toward 
HE 1341-1020 

There are 3 Mg ii systems with Fe ii along this sight hne. 
All these three systems are strong enough to be included 
in our sample. 

5.1.1. Zabs= 0.8728 system toward HE 1341-1020 



= 0.8728 = 1.2778 z, = 1.9154 

abs abs ,— , aba 



3 



Z 




-60 -30 30 -20 20 -40 -20 20 

Relative velocity (kms"') 



Fig. 7. Voigt profile fits to the absorption systems seen 
along the line of sight toward Zcm= 2.135 QSO HE 1341- 
1020. The vertical dotted lines mark the locations of sub- 
components. 

The Voigt profile fits to some of the lines from this 
system are shown in Fig. This system shows Al ii, 
Mg I A2852, Ca iiAA3934,3969 lines in addition to Mg ii 
and Fe ii lines. The Mg ii profile is fitted with 8 compo- 
nents, the three stronger of which show absorption in other 
species. Fe iiA2600 line is blended with absorption lines 
from other intervening systems and the Fe iiA2374 pro- 
file is affected by a bad pixel. In order to avoid unwanted 
systematics we do not consider these lines for the a mea- 
surement analysis. However we use these profiles to check 
the consistency in the number of components and the level 
of saturation in each component. The Al ii line is affected 
by noise and Ca ii lines are too weak to give any meaning- 
ful constraints. We therefore use Fe iiAAA2344,2382,2386, 
Mg iiAA2796,2803 and Mg iA2586. 

5.1.2. Zabs= 1.2778 system toward HE 1341-1020 

The absorption profile of thi-s system is spread over 200 
km s~^. Absorption lines of Mg i, Mg ii, Al ii, Al ill, Fe ii, 
Ca II, Si II, Mn ii, Ni ii, and Ca ii are detected. The core 



is strongly blended over 100 km s^^but a well detached 
narrow satellite is seen ~100 km s~^ away. We have used 
Fe II and Mg ii absorptions from this sharp component for 
the measurement of Aa/a . We do not use Mn ii, Ni ii 
and Ca ii because they are weak. Si ii and Al ii lines fall in 
the Lya forest and the profile of Fe iiA2586 is affected by 
the wing of a contaminating neighboring line. The satellite 
component is well fitted with two components (see mid- 
dle panel in Fig.Cj). Even though the two components are 
merged, we notice that the components are separated by 
more than the velocity dispersion of the individual com- 
ponent. 

5.1.3. Zabs= 1.9154 system toward HE 1341-1020 

This system shows absorption lines from Mg ii, Mg i, Al ii, 
Al III, C IV, Si II, and Si iv. High-ionization lines are 
weak and their profiles are very different from that of low- 
ionization lines. Mg ii lines are saturated and blended with 
atmospheric absorption lines. Thus, we use Si iiA1526, 
Si 11AI8O8 and Fe 11 lines for Aa/a measurement. We 
do not use Fe iiA2586 however as its profile is affected 
by noise. The system is fitted with 4 distinct components 
that are clearly visible in the profiles of the weak lines (see 
right panel in Fig.[7|). Even though the overall fit is good, 
one can notice that the Si 11AI8O8 line profile is apparently 
not reproduced well. This is due to the effect of noise in 
a single pixel in the red wing of the strongest absorption 
component. This shows how delicate the measurement can 
be. 

5.2. Systems along the line of sight toward 
Q 0122- 380 

There are 4 Mg 11 systems along this line of sight. The 
system at 2;abs= 1-9102 has weak Fe 11 absorption lines 
and do not fulfill our criteria to enter our main sample. 

5.2.1. Zabs= 0.8221 system toward Q 0122-380 

This system is defined by narrow Mg 11, Mg i and Fe 11 
absorption lines. All the lines are well fitted by a single 
component. An extra component in the red wing is re- 
quired for the Mg 11 doublets. The Voigt profile fits to the 
system are shown in Fig. |H1 

5.2.2. Zabs= 0.8593 system toward Q 0122-380 

This system is defined by Mg 11, Mg i, Al 11, Al ill, Ca 11, 
and Fe 11 absorption lines. As the spectral range in which 
Al II absorption is redshifted is of lower signal to noise and 
Ca II doublets and Fe iiA2474 are weak we have not used 
these lines for the Aa/a analysis. Resulting Voigt profile 
fits are shown in Fig.|Sl We notice that the Fe iiA2586 line 
is slightly under-predicted in our best model fit. However 
other three Fe 11 lines are very well fitted. We also notice 
that P^r degree of freedom is very good when using 
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Fig. 8. Voigt profile fits to the four absorption systems seen along the line of sight toward Zcm= 2.190 QSO Q 0122—380. 
The vertical dotted lines mark the locations of subcomponents. 



different b parameters for Fe ii and Mg ii. The result on 
Aa/a presented for this system is therefore obtained with 
different b parameters for Mg ii and Fe ii lines. 



5.2.3. Zabs= 1.2433 system toward Q 0122-380 



Mg II absorption profiles in this system is spread over 
~ 70 km s~^. The components are heavily blended. The 
strongest Mg ii component shows Mg i, Al ill, Si ii, Si iv, 
C IV, Ca II, Ni II and Fe ii. Fe iiA2586 and Fe iiA2600 
lines are not covered by our spectrum. Three Fe ii lines 
and the Mg ii doublet can be fitted with 4 components. 
Ni IIA1704 line is blended and Ni iiA1741 is weak and falls 
in the noizy region of the spectrum. We notice small profile 
inconsistencies between the Ca ii doublets. Thus we avoid 
using Ni ii and Ca ii lines in the analysis. The results are 
summarized in Fig. |H1 



5.3. Systems along the line of sight toward PKS 
0237-23 

This QSO is known for the presence of a super cluster of 
C IV absorption lines at Zabs^ 1-6. Our spectrum reveals 5 
Mg II systems one of which at 2;abs= 1.1846 has very weak 
Fe II line and is not considered for our main analysis. 

5.3.1. Zabs= 1-3650 toward PKS 0237-23 

This system shows absorption lines of Fe ii, Al ii, Al III, 
C I, C I*, C II, C II*, Si II, Si IV, Ca ii, Mg i and Mg II 
lines. The Lya line is just below our wavelength coverage. 
The absorption profiles of strong resonance lines spread 
over 230 km s^^. The main component has large num- 
ber of sub-components heavily blended. There is a well 
detached satellite component that could be used for the 
Aa/a measurements. Presence of C i suggests that this 
could possibly be a sub-DLA. In such systems inhomo- 
geneities are generic. Thus it is better to avoid them in 
the Aa/a measurements. We notice that this system hap- 
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Fig. 9. Voigt profile fits to the five absorption systems seen along the line of sight toward Zcm= 2.222 PKS 0237—23. 
The vertical dotted lines mark the locations of subcomponents. 



pens to be in common with our data and that of Murphy 
et al.(2003). Thus we perform the Aa/a measurement in 
the satellite components. We use Al ii, all Fe ii lines (apart 
from Fe iiA2383), and Mg ii lines for our analysis. Murphy 
et al (2003) has found Aa/a — -0.19±0.5 for this system. 
Unlike our case Murphy et al have used the whole pro- 
file of Al II, Si II, Al III, and two Fe ii lines. Our analysis 
gives the best fitted value of O.OiO.l. This is consistent 
with that measured by Murphy et al and our final result. 
The increased accuracy in our measurement is contributed 
by the enhanced S/N and our choice of the well detached 
and unblended satellite component for the analysis. For 
the reason mentioned above we have avoided this system 
in our analysis. 

5.3.2. Zabs= 1-6371 toward PKS 0237-23 

The absorption profiles of this system is spread over ~280 
km s~^. Absorption lines from Mg i, Mg ii, Al ii, Al ill 
Si II and Fe ii are detected. Even though C IV and Si IV 



absorption lines are detected the system is dominated by 
low ionization species. There are two distinct velocity com- 
ponents at Zahs= 1.63588 and 1.63715 that we use for 
Aa/a measurements. The Voigt profile fits to these sys- 
tems are shown in Fig. |51 



5.3.3. Zabs= 1-6574 toward PKS 0237-23 



This system shows absorption lines from Mg ii, Mg i, Al ii, 
Al III, Si II, C IV, Si IV in addition to Fe ii. The low 
and high-ionization lines show distinct profiles. C iv and 
Si IV absorption lines are spread over 400 kni/s. We use 
the uncontaminated Fe ii lines with Si ii and Al ii lines 
as anchors. Mg ii line is heavily blended and has large 
number of extra components compared to Fe ii and Mg i 
is weak. The results of the Voigt profile fits are shown in 
Fig. El 
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5.3.4. Zabs= 1.6724 toward PKS 0237-23 

This is one of the few sub-DLAs (damped Lya systems) 
that are detected in the large programme sample. Mg i, 
Mg II, C i,C I*, C ii,C II*, Si II, Si IV, C iv, Ni ii, Zn ii, 
Cr II, Mn II and Fe ii lines are present. Most of the stan- 
dard resonance lines are saturated. Zn ii, Mn ii, Cr ii 
and Ni ii absorption lines are weak and therefore will not 
give good constraints. The presence of the whole range 
of ionization states and of fine-structure absorption lines 
strongly suggests strong inhomogeneities in the gas. To 
keep our analysis unbiased we did not use this sub-DLA 
system in our analysis. 

5.4. Systems along the line of sight toward 
HE 0001-2300 

There are 5 Mg ii systems detected along this sight line. 
Two of them at Zabs= 0.9489 and 1.6517 are weak systems 
and are not considered in our main sample. 

5.4.1. Zabs= 0.4524 system toward HE 0001-2300 

This system has two very narrow components separated 
by ~ 70 km s~^. The low redshift component shows Mg ii, 
Mg I, Fe II and Ca iiA3934 absorption lines. The high red- 
shift component has very weak Fe ii lines. Thus we use 
only the low redshift component for Aa/a measurement. 
Voigt profile fits to the absorption lines from this compo- 
nent are shown in Fig. ^ Fe iiA2344 and Fe iiA2586 are 
blended with other absorption lines and were not used in 
the analysis. 

5.4.2. Zabs= 1.5855 system toward HE 0001-2300 

Like the 2:abs= 0.4524 system this system has two well 
detached components separated by ~ 110 km s~^. Mg ii, 
Mg I, Fe II Al II, Al III, C iv, Si ii and Si iv absorption 
lines are seen in both components. However Fe ii lines 
are strong only in the high-z component. We notice that 
Mg IIA2803 is contaminated by an atmospheric line. As 
Si IIA1526 is also blended, we use Al II and Mg i absorption 
lines as anchors. Even though the overall fit to the data is 
good, our profiles consistently over predict the blue wing 
of Fe IIA2374. This is probably due to one single pixel in 
the wing being affected by noise. As the two components 
that are required to fit the profiles are very close to each 
other this system is not considered for the final analysis. 

5.4.3. Zabs= 2.1839 system toward HE 0001-2300 

Absorption profiles produced by this system are spread 
over 400 km s~^ Mg i, Mg ii, Al ii, Al ill. Si ii, Si iv and 
C IV absorption lines are detected. Fe ii lines are strong 
only in two of these components (at Zabs^ 2.1854 and 
zabs= 2.1872). As Fe iiA2383, Fe iiA2586 and Fe iiA2600 
absorption lines fall in the wavelength range affected by 



atmospheric absorption, we consider only the clean lines 
in each case. Results of profile fitting are shown in Fig. ^| 

5.5. Systems along the line of sight toward 
QO 109- 3518 

There are two Mg ii system seen along the line of sight 
toward this QSO. We used both of them for Aa/a mea- 
surement. 

5.5.1. Zabs= 1-1827 system toward Q0109-3518 

This system shows absorption lines from Mg ii, Mg i and 
Fe II. The Fe ii and Mg i lines are fitted with a single 
component but Mg ii requires additional components in 
the red wing. Results of Voigt profile fitting are shown in 

Fig. nil 

5.5.2. Zabs= 1.3483 system toward Q0109-3518 

This is a strong Mg ii systems with heavily blended ab- 
sorption profiles spread over 300 km s~^. Absorption lines 
from Mg I, Mg II, Fe ii, Al ii, Al ill. Si ii, Si IV and C IV 
are detected. There are two well detached satellite compo- 
nents. We used the Zabs= 1.3483 satellite for Aa/a mea- 
surement. As Mg II is blended. Si ii, Al ii and Mg i lines are 
used as anchors. For the other satellite the anchor lines are 
all blended. Results of Voigt profile fitting are presented 
in Fig. ini 

5.6. Systems along the line of sight toward Zem= 
2.414 HE 2217-2818 

There are five Mg il systems along the line of sight toward 
this QSO. Three of these systems (at Zabs= 0.7862, 1.6277 
and 1.5556) are weak and are not considered in the main 
sample. The other two systems were used for Aa/a mea- 
surement. 

5.6.1. Zabs= 0.9425 system toward HE 2217-2818 

This system shows absorption lines from Mg ii and Fe ii. 
There are two distinct subsystems that are separated by ~ 
200 km s""'^. Only the high redshift subsystem shows Fe ii 
absorption. Fe ii lines are fitted with four well separated 
components and Mg ii requires an extra component in 
the red wing. We notice that the goodness of the fit in this 
system is decided by how well we model the Mg ii profiles. 
Results of profile fitting arc summarized in Fig. 

5.6.2. Zabs= 1.5556 system toward HE 2217-2818 

This system shows absorption due to Mg ii, Al ii. Si ii, 
Si IV, C IV and Fe ii. The Mg ii absorption profiles are 
spread over ~ 100km s~^. However only the central ve- 
locity component shows other singly ionized species (see 
Fig. Ill|l . Surprisingly this component also shows absorp- 
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Fig. 10. Voigt profile fits to tlie five absorption systems seen along tlie line of sight toward 
HE 0001—2300. The vertical dotted lines mark the locations of subcomponents. 



2.263 QSO 



tion due to O i. High ionization lines such as C IV, Si iv 
are weak and Al ill is absent. 



5.8. Systems along the line of sight toward 
Q 0453-423 



5.7. Systems along the line of sight toward Zem= 
2.435 QSO Q0329-385 

There are two absorption systems detected along this line 
of sight. The Zabs^ 1.43799 system is weak and is not 
considered in our sample. 

5.7.1. Zabs= 0.7627 system toward Q0329-385 

This system shows absorption lines due to Mg i, Mg ii, 
and Fe ii. The Mg ii absorption profile is spread over 200 
km s~^. However the other lines are detected only in the 
central component. The Voigt profile fits to this compo- 
nent is shown in Fig. ^| The profiles are fitted with two 
components and as the velocity separation is smaller than 
the largest b value of the components we do not consider 
this system in our study. 



There are five systems with detected Fe ii and Mg ii ab- 
sorption lines. As most of the Fe ii absorption lines of 
the 2;abs= 0.72617 system are redshifted in the wavelength 
range contaminated by intergalactic Ly-a H i absorptions 
we do not use this system for Aa/a measurement. The ab- 
sorption profiles of the 2;abs= 1.1492 system is spread over 
400km s^^. The Mg ii and Fe ii lines for which we know 
the q values are all saturated. The profiles of Mg i, Ca ii 
and Mn II suggest a very complex blend. Thus we do not 
consider this system for Aa/a measurement. Finally, the 
system at Zabs^ 1.0394 show only Fe iiA2600 absorption 
and is therefore not considered in the analysis. 
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Fig. 11. Voigt profile fits to two of tlie absorption systems seen along the line of sight toward Zem= 2.404 Q0109— 3518 
(first and second column from the left) and Zom= 2.414 HE 2217— 2818 (columns 3 and 4 from the left). The vertical 
dotted lines mark the locations of subcomponents. 



5.9. 2;abs= 0.9083 system along the line of sight 
toward Q 0453- 423 

This system shows detectable absorption from Mg i, Mg 11, 
Fe II, Ca II and Al 11. Fe iiA2586 and Fe iiA2374 are con- 
taminated by absorptions from other systems. Based on 
the Fe 11 profiles the system is decomposed in to eight sub- 
components. The profile fits to the Mg 11 and Fe 11 lines 
free of contamination are shown in Fig. 1121 

5.10. 2abs= 1.8583 system along the line of sight 
toward Q 0453- 423 

This is a narrow single component absorption system de- 
tected in Mg II, Al II, Al III, Si 11, Si iv and Fe ii. Results of 
the single component Voigt profile fit to three prominent 
Fe II lines and the Mg 11 doublet are shown in Fig. El 



5.11. Zabs= 2.3004 system along the line of sight 
toward Q 0453- 423 



This system shows absorption lines of Mg 11, Fe 11 and 
Al II. Fe IIA2586 and Fe iiA2600 are not covered by our 
spectrum. Fe iiA2374 is very weak. The signal-to-noise ra- 
tio in the spectral range where other Fe 11 lines are seen is 
^ 30 only. In addition, both Mg 11 lines are contaminated 
by atmospheric absorption. There is also profile inconsis- 
tencies between the weak Al 11 lines (a potential anchor) 
and Fe 11 lines. Thus we do not consider this system in our 
study. 
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Fig. 12. Voigt profile fits to tlie absorption systems seen along the line of sight toward Zem= 2.435 QSO Q0329— 385 
(first column from the left), Zem= 2.611 QSO HE 1347 - 2457 (column 2 from left), and Zem= 2.658 Q 0453 - 423 
(columns 3 and 4 from the left) . The vertical dotted lines mark the locations of subcomponents. 



5.12. Absorption systems along the line of sight 
toward Zcm= 2.767 PKS 0002- 422 

There are 5 Mg ii systems with detected Fe ii absorption 
lines along this line of sight. The Zabs= 0.8367 system is 
a very strong Mg ii system whose absorption profile is 
spread over '-^600 km s~^. As the corresponding absorp- 
tion profiles are heavily blended we do not consider this 
system in our analysis as is the case for the Zabs= 1.9885 
system in which Fe ii lines are too weak. The remaining 3 
systems are considered in our main analysis. 

5.12.1. Zabs= 1.5418 system toward Q 0002-422 

This system is defined by Mg ii, Mg i, Al ii, Al ill, and Fe ii 
absorption lines. The Mg ii absorption lines are spread 
over ~180 km s~^. The other low ionization lines are seen 



only in the strongest central component. Voigt profile fits 
to the absorption lines are shown in Fig. ^] 

5.12.2. Zabs= 2.1678 system toward Q 0002-422 

This system is defined by Mg ii, Mg i, Al ii, Al ill, Fe ii, 
Si II, Si IV and C iv absorption lines. Like most systems 
in our sample the high and low ionization profiles have 
different profiles, the Mg ii absorption profile is spread 
over 80 km s^^. However only the strongest central com- 
ponent shows other low ionization absorption lines. As the 
central Mg ii component is blended in the red wing with 
other components we use Si ii and Al ii as anchors. The 
presence of Fe 11AI6O8 with negative q coefficient makes 
this system very sensitive to Aa/a . The system is fitted 
with two components and the velocity separation between 
them is larger then individual b values. 
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Fig. 13. Voigt profile fits to two of tlie absorption systems seen along the line of sight toward Zcm= 2.767 Q 0002—422 
(First three columns from the left) and Zom= 3.280 QSO PKS 2126—158 (last column). Vertical dotted lines mark the 
locations of subcomponents. 



5.12.3. Zabs= 2.3008 system toward Q 0002-422 



This system is defined by Mg ii, Mg i, Al ii, Al ill, Fe ii, 
Si II, Si IV and C iv absorption lines. Like most systems 
in our sample the high and low ionization profiles have 
different profiles, the Mg ii absorption profile is spread 
over 200 km s^^. The main component is complex and 
blended. There is a well detached satellite in the blue wing 
of the profile. We use this for our Aa/a measurement. We 
use the unblended Fc ii lines and Si iiA1526 and Al ii lines 
for our analysis. The three components that fit the profile 
are well resolved. 



5.13. Absorption systems along the line of sight 
toward Zem= 3.280 QSO PKS 2126-158 

Two Mg II systems are detected along the line of sight 
toward this QSO. The Zabs= 0.6631 system shows absorp- 
tion fines from Mg II, Mg I, Fe Ii, and Ca Ii. AU the fines 
are heavily blended with intervening Lya absorption lines 
owing to the large redshift of the QSO. We do not consider 
this system for our analysis. 

5.13.1. Zabs= 2.0225 system toward PKS 2126-158 

This system shows absorption due to Mg ii, Si ii and Fe ii. 
The Mg II profile is spread over 140 km s~^. The strongest 
Mg II component is saturated. We use the narrow satellite 
component in the low redshift side, at Zabs= 2.02192 for 
Aa/a measurement. As the Mg iiA2796 is affected in the 
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Table 4. Results for the weak systems that show at least 
two detectable Fe ii lines: 



QSO 




Aa/a (in 10"^) 


Q 0122 - 380 


1.9102 


-0.9 ± 1.6 


PKS 1448 - 232 


1.5847 


+0.4 ±2.0 


PKS 0237 - 23 


1.1846 


-0.9 ±2.7 


HE 0001 - 234 


0.9489 


-1.1 ± 2.3 


HE 2217- 2818 


1.6277 


+0.3 ±0.6 


Q 0329 - 385 


1.4379 


-0.4 ±0.6 


HE 1347 - 2457 


1.5082 


-3.9 ± 1.4 


PKS 0329 - 255 


0.9926 


-0.3 ± 1.3 


HE 2347 - 4342 


1.7962 


-0.4 ±2.0 



blue wing, we use Si II and Mg ilA2803 as anchors. Voigt 
profile fits to the lines used in the analysis are shown in 

Fig.m 

6. Results 

6.1. Fitting procedure 

As discussed before the accuracy of the method depends 
on how well absorption line profiles are modeled using 
Voigt profiles. As can be seen from Table. ^ the relative 
shifts between different Fe ii lines (apart from Fe iiA1607) 
are small compared to the shifts between the Fe ii lines 
and the anchors. We use this fact to refine our fitting pro- 
cedure. 

We first fit all the Fe ii lines simultaneously using labo- 
ratory wavelengths {Aa/a = 0). This allows us to find out 
about (i) bad pixels, (ii) unknown contaminations and (iii) 
the velocity component structure in Fe ii. We therefore can 
remove inadequate profiles (because of contamination or 
bad pixels) from the Aa/a analysis. Similar exercise was 
carried out for Mg ii doublets and other anchors. Based on 
these preliminary fits a first set of parameters is generated 
to start the Voigt profile fitting procedure that includes 
Aa/a variations as described in section 3.3. 

In all our fits we have tied up the redshifts of all species 
under consideration. We do not know a priori whether 
b parameters of different species are identical or not. To 
take this uncertainty into account we perform the fit of 
the systems in two cases: (i) assuming that b parameters 
of individual components are same for all species and (ii) 
assuming different b parameters for different species. We 
determine the best fitted Aa/a values for both cases. It is 
found that in all cases the two values are consistent with 
one another within la uncertainty. For our final result we 
consider the fit that gives the smaller reduced x^. 



6.2. Determination of Aa/a 

The analysis for individual systems is shown in Fig. 1141 
Here, we plot as a function of Aa/a . The dot with 
error-bar in each panel gives our best fit estimate of 
Aa/a and its errors. The vertical lines in each panel shows 
the value Aa/a = 0. It is clear from the figure that apart 



from few cases the curve is well behaved and the best fitted 
values are consistent with zero in most cases. The numer- 
ical values of the fits are given in Table. 01 In this table 
we list in the first four columns the name of the QSO, the 
emission redshift (zem), absorption redshifts (zabs) and ve- 
locity dispersions b of individual components. When the 
best fit solutions are obtained with different b parameters 
the b values quoted are for the main anchor line (Mg ii. 
Si II or Al II lines). Column #5 lists the transitions used 
in the analysis with the same notations as in Tabled A''c 
given in column #6 is the number of subcomponents used 
to fit the Fe ii absorption profiles. The last two columns 
give the determination of Aa/a obtained using labora- 
tory wavelengths of Mg ii, Mg I, and Si II absorption lines 
with terrestrial isotopic abundances (case 1) and with the 
wavelengths of the dominant isotopes (case 2, see next 
Section). Numbers in brackets are the reduced for the 
fit. The numbers marked with a asterisks (*) are values 
obtained using different b values for different species. 

In Fig. El we plot individual determinations of 
Aa/a as a function of Zabs- We also plot the exisiting re- 
sults from the literature. The horizontal dotted line gives 
the weighted mean of our sample with 1/error^ weight- 
ing. We estimate the error in the weighted mean using the 
standard equation. 



Error in x„ 



(N - l)Sf Wi 



(2) 



Here, is the weighted mean of the variable xi of 
sample size N whose weights are Wj. The shaded region 
passing through most of the error bars is our measured 
weighted mean and its 3a error. The histogram in the left 
hand side of the panel shows the distribution of Aa/a. 
The weighted mean with 1 /cr^ weighting obtained for our 
sample is — 0.06±0.06x 10~^ and the standard deviation in 
our measurements around the mean is 0.41x10^'''. All the 
points used are consistent with this weighed mean value 
with a reduced oi 0.95. As can be seen from the figure 
there are two points that deviate by more than la from 
derived weighted mean value. These are from the 2:abs= 
0.9425 system toward HE 2217-2818 and the Zabs= 1.9154 
system toward HE 1341—1020 (see also the curve for 
these systems in Fig. ll4|l . As can be seen from Fig.^] the 
Fe II lines in the 2;abs= 0.9425 system are weak (barely 
above the cutoff). The reduced x^ is dominated by un- 
certainties in fitting the Mg II profiles. In the case of the 
latter system the effect could just be due to the odd pixel 
in the Si 11AI8O8 close to u = km s"^(see Fig.jTj). 

Just for completeness we also fitted the weak systems 
that show at least two Fe 11 lines. The recovered values 
of Aa/a based on these are given in Table 01 As ex- 
pected the individual measurements have large errors. The 
weighted mean value of Aa/a measured from these sys- 
tems is — 0.40±0.36. Most of the individual measurements 
and the weighted mean are consistent with the value ob- 
tained from our main sample. 
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Table 3. Rcxiovc^rcxl value' of Ao/fi in individual systc'ms 



Name 








^abs 






b 




Transitions 


No 


Ack/ck (in un 
case 1 


its of 10"-'') 
case 2 


HE1341 


-1020 


2.135 


0, 




_j_ Q 


,000046 


5, 


.65 it 0, 


13 


abdfgi 


3 


+0.0 ± 0.2(1.19) 


-0.6 ± 0.4(1.20) 











.872841 


± 


.000212 


5 


.22 ± 0, 


,28 


abdfg 















.cs / ^yuo 


i 


,000195 


3, 


.28 ± 0, 


,28 


abdfg 



















i 


,000266 


2 


.49 ± 0. 


,90 


ab 

















S79499 


-j- Q 


,000089 


2 


.59 ± 0, 


,30 


ab 














0, 


.872477 


± 0. 


,000563 




2.68 




ab 














0, 


R79^RR 


± 


,000354 


7, 


.65 ± 0, 


,84 


ab 

















. o / OUU / 


± 


,000206 




0.50 




ab 














1 


277834 


it 


,000584 


4, 


.03 it 0, 


,51 


abdefh 


2 


+0,1 ± 0.2(1.01) 


+0.0 ± 0.50(0.97) 








1 


.277780 


± 0. 


,000571 


3 


.35 ± 0, 


,52 


abdefh 














1 


915216 


rh 


,000628 


2 


.57 ± 0, 


,73 


cdefhkl 


4 


+0.8 ± 0.3(1.49) 


+0.6 ± 0.3(1.48) 








1 




db 0. 


.000226 


3 


.40 ± 0. 


.23 


cdefhkl 














1 


Q1 ^4R^ 


± 


,000302 


6, 


.24 ± 0, 


.28 


cdefhkl 














1 


.915561 


± 


002269 


3 


.86 ± 1 


,44 


cdefhkl 








Q0122- 


380 


2.190 





R99n7fi 


i 


000032 


5 


.69 ± 


,07 


abdfghi 


1 


+0.0 ± 0.9(0.87) 


+0.0 ± 0.9(0.87) 











S991 S9 


dr 


,000457 


6, 


.41 ± 1 


,28 


ab 



















i 


.000211 


3, 


.92 ± 0. 


.08 


abfghi 


2 


-0.3 ± 0.2(1.29) 


-0.6 ± 0.3(1.18) 











.csoyoiis 


_j_ Q 


,000646 


3, 


.42 ± 0, 


,13 


abfghi 














1 


.243293 


it 0. 


,000109 


5 


.09 ± 0, 


,12 


abdef 


3 


-0.1 ± 0.1(0.89) 


-0.6 ± 0.1(1.01) 








1 




i 


000344 


3, 


,12 ± 


50 


abdcf 














1 


243483 


i 


000122 


6, 


,83 ± 0, 


41 


abdcf 














1 


243589 


i 


,000731 


5, 


,20 ± 1 


,03 


ab 








PKS0237-23 


2.222 


1 


fi'^^f;74 


it 


,000708 


9, 


.23 it 1, 


,06 


abk 


2 


+0.2 ± 0.7(0.82)* 


-0.2 ± 0.3(0.86)* 








1 


R^=;7Qn 


-j- Q 


,000132 


2 


.56 ± 0. 


,34 


abdghk 














1, 


. OOOOoO 


i 


.000050 


4, 


.04 ± 0. 


.80 


abdghk 














1 


. DODU'ii- 


-j- Q 


,000152 


4, 


.20 ± 0, 


,33 


abk 














1 


R^71 9R 


it 


000114 


3 


,60 it 


42 


abcdgh 


6 


+0.6 ± 0.6(1.16)* 


-0.4 ± 1.3(1.19)* 








1 


R3fiQn3 


it 0. 


000119 





,40 ± 0, 


13 


abdgh 














1 


.DO / OOi- 


-j- Q 


,000244 


3, 


.87 ± 0, 


,62 


abdgh 














1 


.DO / J-OO 


i 


,000209 


15.4 ± 0. 


.41 


abcdgh 














1 


.637437 


± 0. 


,000170 


6, 


.17 ± 0, 


,24 


abdgh 














1 


657390 


it 


,000223 


5.46 it 0, 


,34 


dfjk 


8 


+0.3 ± 0.5(0.92) 


+0.3 ± 0.7(0.92) 








1 




i 


,000932 


5, 


.51 ± 0, 


,92 


dfjk 














1 


657415 


it 


,000348 


11.57 ± 0.34 


dfjk 














1 


R^RRR.1 
.DODOD^ 


-j- Q 


.001122 


10.42 ± 1.22 


fik 














1 


RK;7nRR 
.DO ( UDD 


_j_ Q 


,001358 


4, 


.75 ± 1. 


,34 


fjk 














1 


R^779^ 
.DO / / ZO 


± 


,000785 


6, 


,91 it 1, 


,02 


tjk 














1 


R^7SR4 
.DO / Oo'i 


± 


,000262 




,36 ± 0, 


28 


dfjk 














1 


.DO I oy 1 


it 


,000622 


2 


.16 it 0, 


,17 


dfjk 










2340 


2.263 





zLf^9nfin 

.'iOZlJDD 


i 


000024 


1 


.74 it 0, 


,05 


abefhim 


1 


+0.2 ± 0.5(1.10) 


-0.3 ± 0.6(1.21) 








2 




i 


001503 


2 


.03 it 0, 


,23 


bdegjk 


3 


-0.2 ± 0.2(1.20) 


-0.2 ± 0.3(1.31) 








2 


187237 


-j- Q 


,000129 


4, 


.73 ± 1. 


,06 
















2, 


187106 


-j- Q 


.001420 


6 


.61 ± 0. 


.66 
















2 


185294 


-j- Q 


,000089 


2 


.29 ± 0. 


,19 


abdefjk 


2 


+0.2 ± 0.3(1.15) 


-0.2 ± 0.3(1.17) 








2 


.185452 


it 0. 


,000264 


3, 


.48 ± 0, 


,52 


abdefjk 














2 


. lOOZOD 


it 


000512 




1.51 




abjk 














2 


. i-OOOOZ 


-j- Q 


,000839 




1.77 




abjk 










3518 


2.404 


1 


1 R9RR^ 


i 


.000039 


3, 


.04 ± 0. 


.09 


abdfhi 


1 


+0.0 ± 0.8(0.98) 


-0.6 ± 0.6(0.98) 








1 


182770 


± 


,000389 


3, 


.17 ± 1, 


,18 


ab 














1 


.348308 


it 0, 


,000204 


2 


,28 it 0, 


,21 


defghijk 


2 


-0.6 ± 0.4(1.08) 


-0.6 ± 0.4(1.08) 








1 


.348346 


± 0, 


000415 




3.72 




dcfghijk 








HE2217 


-2818 


2.414 





.942330 


± 0, 


000058 


4 


,91 ± 


12 


abfh 


5 


-1.2 ± 0.7(0.90)* 


-1.6 ± 0.8(1.36) 











.942574 


± 0, 


()0()067 


4 


,79 ± 


15 


abfh 

















.942655 


± 0, 


000119 


3 


,27 ± 0, 


75 


abfh 

















.942701 


± 0. 


,002092 


8, 


.88 ± 2. 


,20 


abfh 














0, 


.942435 


±0. 


.000464 


7, 


.38 ± 0. 


,64 


abfh 














0, 


.942476 


±0, 


,000282 




4.82 




abfh 














1 


.555425 


it 0, 


,000029 


2 


,52 it 0, 


,06 


ab 


1 


+0.2 ± 0.6(1.22)* 


+0.0 ± 0.3(1.35)* 








1 


.555580 


it 0, 


,000098 


2 


.35 ± 0, 


,06 


ab 














1 


.555721 


± 0, 


,000178 


3, 


.08 ± 0, 


,39 


ab 














1 


.555886 


±0. 


.000041 


4, 


.98 ± 0. 


,07 


abdfghk 














1 


.556162 


±0, 


,000055 


3, 


.76 ± 0, 


,11 


ab 








HB1347 


-2457 


2.611 




1.439541 




40.91 




k 


3 


-0.0 ± 0.6(1.10) 


-0.4 ± 0.6(1.11) 








1 


.439219 ± 0.002721 


14.2 ± 


90 


dgkl 














1 


.439243 ± 0.000365 


6, 


.96 it 0, 


,36 


dgkl 














1 


.439349 ± 0.000195 


2.61 ± 0. 


,07 


dgkl 
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Fig. 14. Results from the UVES sample: is plotted as a function of Aa/a for each system in our sample. The 
minima of the curve (marked with a dot) gives the best fitted value Aa/a . The error in the measurement (error bar 
around the dot) is derived using Ax^ = 1- Name of QSOs and redshifts Zabsare given in each panel. The vertical dotted 
line in each panel marks the location of Aa/a — 0. It is interesting to note that most of the systems are consistent 
with Aa/a = within error bars. Two of the systems Zabs^ 1.5864 toward HEOOOl— 2340 and Zabs= 0.7627 toward 
Q 0329—385 are not considered in the analysis as the only two components in these systems are blended as per our 
definision. 



A summary of the results for different sub-samples are 
given in Table The sample identification is given in the 
first column. Number of systems used and the median red- 
shift of the sample are given in column #2 and #3. The 
mean, weighted mean and cr of the measured Aa/a val- 
ues are given in columns #4, #5 and #6 respectively. Last 
column gives the reduced Xw obtained from all our mea- 
surements for the measured weighted mean. It can be seen 
that in none of our subsamples we find a significant change 
in a. 

6.3. The magnesium isotopic abundance 

One major uncertainty in the many-multiplet analysis 
comes from the determination of the effective rest wave- 



lengths. Even though laboratory wavelengths are mea- 
sured with a precision of a few 0.1 mA, values given 
in Table. Q assume terrestrial isotopic abundances. This 
assumption may not be valid at high redshift. In the 
Astrophysical settings the effect of isotopic shifts could 
be important for Si and Mg (see Murphy et al. 2003, 
Ashenfelter et al. 2003). Indeed, the relative abundances 
of different isotopes may depend on the overall metallic- 
ity of the gas. Gay & Lambert (2000) have shown that 
the abundance of ^^Mg and ^^Mg relative to ^^Mg de- 
creases with decreasing metallicity. In the low metallicity 
gas (Z < O.OIZq) most of the metals will be in their dom- 
inant isotopic state. Thus in the extreme case of very low 
metallicity, the effective rest wavelengths could take val- 
ues in the range from terrestrial composite wavelength to 
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Fig. 15. Results from the UVES sample: The measured values of Aa/a from our sample (filled circles) are plotted 
against the absorption redshifts of Mg ii systems. Each point is the best fitted value obtained for individual systems us- 
ing minimization as demonstrated in Fig. 1141 The open circle and stars are the measurement from Oklo phenomenon 
and from molecular lines respectively. The weighted mean and la range measured by Murphy et al.(2003) are shown 
with the horizontal long dashed lines. Clearly most of our measurements are inconsistent with this range. The shadow 
region marks the weighted mean and its 3ct error obtained from our study [< Aa/a >w = (— 0.06 ± 0.06) x 10^^]. 
Our data gives a 3a constraint on the variation of Aa/a to be —2.5 x 10~^® yr^^ < (Aa/aAt) < 1.2 x 10^^^ yr^^ in 
the case of flat universe with fi^ = 0.7, flm. — 0.3 and Ho = 68 km s~^Mpc^^ for the median rcdshift of 1.55. 



wavelengths corresponding to the dominant isotope. This 
range is less than 0.5 mA for Mg i and Si ii absorption 
lines but is of the order of 1 mA for Mg ii lines. In or- 
der to accommodate this uncertainty we fit the systems 
using the wavelengths of the species from the dominant 
isotope. The measurements are given in the last column 
of Table 3. As expected, using these abundances leads to 
a lower a determination (see Table 5). Note however that 
even in this extreme case, the variation stays smaller than 
what has been claimed from previous studies. Note that 
the assumption of very low metallicity is extreme as: (i) 
the systems in the z range 0.4—2 are more likely to have 
metallicity larger than 0.1 Zq (e.g. Ledoux et al. 2002); (ii) 
the measured mean ratios of Mg^'*:Mg^^:Mg^^ in the cool 



dwarf with metallicity Z = -1.5 to -1.0 Zq is 80:10:10 
(from Table. 1 of Yong et al. 2003). This gives the weighted 
mean wavelengths close to the terrestrial wavelengths, and 
(iii) We also notice that the minimum for the fit in most 
cases (apart from 4 cases) are better when we use the lab- 
oratory wavelengths in our analysis. Therefore, although 
some additional uncertainty and scatter could come from 
the isotopic abundances being different from that of ter- 
restrial composition, our result using the laboratory wave- 
lengths is most probably robust. 

7. Conclusion 

We have applied the MM method to a homogeneous sam- 
ple of 50 Mg II systems observed along 18 QSO lines of 
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Table 5. Summary of results for various sub-samples: 



Sample 




Number of 
systems 


z 


mean 


Aa/a 

weighted mean 


a 


Xw 


Single-)-double(case 


1) 


12 


1.54 


-f0.01±0.15 


-0.08 ±0.07 


0.27 


0.55 


Full sample (case 1) 




23 


1.54 


-0.02 ±0.10 


-0.06 ±0.06 


0.41 


0.95 


Weak(case 1) 




9 


1.51 


-0.80 ±0.58 


-0.40 ±0.36 


1.27 


1.00 


Weak-|-All(case 1) 




32 


1.51 


-0.24 ±0.18 


-0.07 ±0.06 


0.81 


0.96 


Single-|-double(case 


2) 


12 


1.54 


-0.31 ±0.14 


-0.33 ±0.09 


0.33 


0.63 


Full Sample(case 2) 




23 


1.54 


-0.33 ±0.11 


-0.36 ±0.06 


0.44 


1.06 



Case 1: Laboratory wavelengths given in Table. □ are used. 



Case 2: Rest wavelengths of dominant isotopes for Mg and Si are used. 



sight observed with UVES at the VLT. Using extensive 
simulations we show our Voigt profile fitting procedure 
works well to recover the input value of Aa/a for sim- 
ple single component systems. We can recover the input 
Aa/a values with an error of fi = 0.23 x 10~^ in the case 
of individual components for data of quality comparable 
to the UVES data (S/N 70). We show that the uncer- 
tainty is about twice larger for data of S/N ~ 30 as used by 
previous surveys. In addition at these S/N ratio, the use 
of strongly blended systems leads to larger uncertainties 
as well. We show that weak lines should be avoided in the 
determination as their use can lead to false determination 
of Aa/a . 

We device the selection criteria for our sample using 
extensive simulations. This leads us to 23 systems avoiding 
stongly blended systems and systems with weak absorp- 
tion lines. Individual systems and fits to all absorption 
lines are discussed in detail. We also avoided the two sub- 
DLAs from the analysis. In one of the systems that is in 
common with Murphy et al. (2003) our results are con- 
sistent with the earlier measurement albeit with smaller 
error (see section 5.3.1). 

The weighted mean value of the variation 
in the fine-structure constant obtained from our 
analysis over the redshift range 0.4<z<2.3 is 
Aa/a -0.06 ± 0.06 X 10"^ For the median 

redshift of our sample (z = 1.54) we obtain a 
3cr constraint on the variation of Aa/a to be 
-2.5 X 10-1^ yr-i < (Aa/aAt) < 1.2 x lO'^^ yr'^ in the 
case of flat universe with Q.x — 0.7, SI™ = 0.3 and Hp = 
68 km s^^Mpc^^. To our knowledge this is the strongest 
constraint from QSO absorption line studies till date. 

We show the effect of varying the isotopic abundances. 
Only if we assume isotopic abundances for very low metal- 
licity would a show some variation. This is non accept- 
able for (i) assuming that all systems in the redshift range 
0.4—2 have extremely low metallicity is incorrect as met- 
als are seen in these systems and metallicity is measured 
to be larger than 0.1 Tjq (e.g. Ledoux et al. 2002); (ii) at 
Z ^ 0.1 Z0, which would be a more reasonable assump- 
tion, the composite wavelength we get will be close to the 
laboratory value we use in our study (we use Yong et al. 
2003 for isotopic abundances). Therefore, although some 
additional uncertainty and scatter could come from the 
isotopic abundances being different from terrestrial, our 



result using laboratory wavelengths is most probably ro- 
bust. 

In conclusion, our study does not support claims by 
previous authors of a statistically significant change in 
Aa/a with cosmic time at z > 0.5. Our result still does 
allow smaller variations in excess of what is found based 
on the Oklo phenomenon. Future very high resolution 
(R~100,000) spectroscopic studies are needed to probe 
the variations in a with much better accuracy. 
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